Charm production from photo-nucleon reaction in a hadronic model 
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We study the total cross section for photo production of charmed hadrons near threshold using 
a hadronic Lagrangian. Both two-body final states involving Ac and a charmed meson as well as 
three-body final states involving nucleon and a charm-anticharm meson pair are included. With 
appropriate cut-off parameters in the form factors at interaction vertices, the model gives a total 
cross section that is consistent with the measured data at center-of-mass energy of 6 GeV. The result 
is compared with the prediction from the leading-order perturbative QCD. 
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I. INTRODUCTION 



Reliable estimates of the production and scattering 
cross sections of open and/or hidden charmed hadrons 
in hadronic matter are important for understanding 
many phenomena in relativistic heavy ion coUisions and 
hadron-nucleus reactions. In particular, the dissociation 
cross sections of J/4' by light meson or nucleon are di- 
rectly related to the role of hadronic suppression of J/^p 
production in heavy ion collisions [1,2], and have thus 
been a subject under active investigations. In the heavy 
charm quark mass limit, a formula based on the leading- 
order (LO) perturbative QCD was derived in Refs. [3,4] 
for the absorption cross sections of J/ip by hadrons. 
Higher-order corrections due to the target mass and rel- 
ativistic effects have also been estimated [5-7]. However, 
the perturbative QCD approach is useful only at very 
high energy when nontrivial higher twist and higher as 
corrections are small. At low energies, these corrections 
become large, so nonperturbative approaches are needed. 
Various phenomenological models have thus been intro- 
duced for studying the J/4' absorption cross section by 
hadrons at low energies. These include the hadronic 
model based on effective Lagrangians [8-10], the QCD 
sum rules [11], and the quark-exchange model [12]. Re- 
sults from these phenomenological models all give much 
larger cross sections for J/tp absorption by light mesons 
than that given by the perturbative QCD formula. To 
test the prediction from the hadronic model, the same 
effective Lagrangian has been used to evaluate the ab- 
sorption cross section of J/tp by nucleon, and it is found 
to be consistent with that extracted from photonuclear 
production of J/ip [9]. 

In the present paper, we generalize the effective 
hadronic Lagrangian, which has also been used to 
study charmed meson scattering by hadrons [13,14] and 
charmed meson production from meson-nucleon scatter- 
ing [15], to include the photon and to study charmed 
hadron production from photo-nucleon reaction near 
threshold. Both two-body {AcD) and three-body (NDD) 
final states are included. We find that using reason- 
able values for the cutoff parameters at the interaction 
vertices, the resulting charmed hadron production cross 



section near threshold energy is consistent with the mea- 
sured one at center-of-mass energy of 6 GeV [16], al- 
though at high energies its value is much smaller than 
that measured experimentally or given by the LO per- 
turbative QCD. We further find that the relative con- 
tribution of two-body to three-body final states is also 
consistent with the experimental data. Our results thus 
provide another independent test and confirmation of the 
validity of hadronic model for determining the cross sec- 
tions involving heavy flavors at low energies. The ef- 
fective hadronic Lagrangian used in the present study 
will also be useful for evaluating other low energy cross 
sections involving heavy flavored hadrons, which can be 
studied at both the Japanese Hadron facility [17] and at 
the GSI future accelerator [18], where open and hidden 
charmed hadrons are copiously produced in proton- and 
antiproton- nucleus reactions near threshold. 

This paper is organized as follows. In Section II, we 
first calculate the cross sections for photoproduction of 
charmed hadrons from nucleon with three particles in the 
final state, which is dominated by pion and rho meson ex- 
change. Photoproduction of charmed hadrons from nu- 
cleon with two particle in the final state, which includes 
the charmed meson exchange, is studied in Section HI. 
In Section IV, we show the total cross section and com- 
pare it to available experimental data and results from 
the LO QCD. Finally, summary and discussions are given 
in Section V. 



II. PHOTOPRODUCTION OF CHARMED 
HADRONS FROM NUCLEON WITH 
THREE-PARTICLE FINAL STATES 



In photoproduction of charmed hadrons from nucleon 
with three particle in the final states, the three possible 
reactions are 77V DDN, jN DD*N{D*DN) ,and 
D*D*N. The lowest-order diagrams for the pro- 
cess jN DDN are shown in Fig. 1 and involves the 
exchange of pion in the intermediate state; those for the 
processes 7iV DD*N{D*DN) and 7iV D*D*N in- 
volve the exchange of rho meson, and the lowest-order 
diagrams for the two processes are shown in 2 and 3. 
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FIG. 1. Photoproduction of charmed hadrons from nucleon 
involving pion exchange. 
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FIG. 2. Photoproduction of charmed hadrons {DD) from 
proton involving rho meson exchange. 



To evaluate the cross sections for these processes, we 
use the same Lagrangian introduced in Refs. [9,14,15] for 
studying charmed meson scattering by hadrons. This La- 
grangian is based on the gauged SU(4) flavor symmetry 
but with empirical masses. The coupling constants arc 
taken, if possible, from empirical information. Other- 
wise, the SU(4) relations are used to relate the unknown 
coupling constants to the known ones. Photon is then in- 
troduced in the Lagrangian via gauging its Uem(l) part 
as in ref. [19]. 
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FIG. 3. Photoproduction of charmed hadrons {D*D*) from 
proton involving rho meson exchange. 

The interaction Lagrangian densities that are relevant 
to the processes shown in Figs. 1, 2, and 3 are given as 
follows: 



JCwNN = -ig-KNNN^^fN ■ TT 
^pNN 



jO-nDD' = i9nDD'D*'^T ■ {Ddf^TT - d^Dlf) + H.C., 

jCpD'D' = igpD'D' [{d^D*-fDl - D*-fd^Dl) ■ ^ 
+ {D*-f-d^fr-d^D*''f-p,)D*>^ 



t^p-^DD 
l^p-yD*D* 



= ieA^[DQd^,D - {d^D)QD], 

= ie[A^'{d^D*''QDl - D^^Qd^Dl) 

+ {d^A-Dl - A''d^Dl)QD*^' 

+ D*^^Q{A^'^^Dl-^^,A-Dl], 

= -eg^DD'A>'{D*^{2TQ - Qt)D 

+ D{2QT-TQ)D;)-n, 

= egpooAi'DiTQ + Qt)D ■ p^, 

= egpD'D'{A''Dt{2TQ - Qr)D; 

+ A-'D;{2rQ~QT)D: 

-A^D*-'{2tQ-Qt)D:)-p^. 



(1) 



In the above, t arc Pauli spin matrices, and tt and p de- 
note the pion and rho meson isospin triplet, respectively, 
while D = {D^,D+) and D* = {D*°,D* + ) denote the 
pseudoscalar and vector charm meson doublets, respec- 
tively. The operator Q is the diagonal charge operator 
with diagonal elements equal to and -1. 

For coupling constants in the above interaction La- 
grangians, we use the following empirical values: g^NN = 
13.5 [20], gpNN = 3.25, Hp = 6.1 [21], g^oD- = 5.56 
[15,22], and 5^z3Z3 = 2.52 [14]. 

The diagrams in Figs. 1,2, and 3 can be separated to 
two types; one in which the photon is coupled to mesons 
such as the first three diagrams (denoted by (ia) to (ic) 
with i=l to 8), and the other in which the photon is 
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coupled directly to the incoming or outgoing charged nu- 
cleon. As shown later, contributions from the latter type 
are much smaller than those from the first type of di- 
agrams and are thus neglected in the following calcula- 
tions. We note that diagrams of the first type are similar 
to those for J/ip absorption on nucleon, which can be 
interpreted as absorption by the virtual pion and rho 
meson cloud of a nucleon. Here, they can be considered 
as charmed hadron production from the nucleon's virtual 
meson cloud. 

The amplitudes for the four processes in Fig. 1 are 
given by 



Mi = -iag^NNN{p3)-/5N{pi)-^ 

X {M,a+Mib + M,^), 



1 



(2) 



with i = 1 to 4, while the amplitudes for the four pro- 
cesses in Fig. 2, 3 can be written as 



Mj = agpNNN(j>3) 



ipi -P3)a 



1 



-9i^u + 



2m N 

(Pi -Pz)ti{Pi-Pz)i^ 



t-mj 



(3) 



with j = 5 to 8. In the above, pi and are four mo- 
menta of the initial and final nucleons, respectively. The 
coefficient a is, respectively, 1 and \/2 for neutral and 
charged pion coupling to nucleon. 

The three amplitudes Mia, Mn,, and M^c are for the 
subprocess ivf — > D*D in Fig. 1. It can be shown that 
they fulfill the chiral constraint [10], i.e., their sum van- 
ishes at the soft pion limit. Explicitly, they are given 

by 

Mia = M2a 

= eg^DD-i-^h + ka)'' 



m 



D 



X {ki - k3 + fc4)''£3^e2,y, 

Mib = M2b 

= -eg^DD'{-ki - k4y 



u — m, 



9af) 



(fci - k4)a{ki - k4)f: 



X [(-fc2 - ksfg'^'^ + (-fci + fc2 + k^Yg^^ 
+ {ki + ks- kiYg^'']eii,e2u, 

Mlc = M2c 

= egTrDD'g^'^£3ti£2i^, 

Msa = V2eg^DD'i-2ki + ksY 



1 



t — mj-i 



X {ki -k3 + kiYe3fj,e2v, 
M3b = -V2eg^DD' (2A;i + fcs)"- 



X (fci + k2 + kiYe3ij,£2u, 
Msc = '2'V2egT,DD*g'^"£s^e2v, 
Mia = -V2egT,DD'{-ki - k^Y 



u — mj-,, 
{ki - k4)a{ki - ki)i3' 



ml. 



X [{-k2 - ksYg'^" + i-ki + k2 + kiYg^"" 

+ (fci-hfc3-fc4)V]^3/.e2>., 

■kdd* (2fci -|- k2Y 



s-mt 



X (fci + fc2 + k4,Y£3^£2u, 

Mac = -V2eg^DD*g'"'£3^i£2i^, 



(4) 



where ki denotes the momentum of particle i of each 
subprocess, and and Si, are the polarization vector of 
7 and D* , respectively. We choose the convention that 
particles 1 and 2 represent initial-state particles while 
particles 3 and 4 represent final-state ones on the left 
and right sides of the diagrams. 

The amplitudes MJ^, Mj^, and Mj^ are those for the 
subprocesses pj DD and pj D*D* in Figs. 2 and 
3, and they are are written explicitly as: 



M^ 



5a 



= -egpooiki - 2/13)'' 

X (fci - fcs + fc4)''e2i., 
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t — m 



D 



■^56 = -egpDD{-ki + 2k4Y 



u — m 



D 



X (-fci - fc3 -I- k4)"£2v, 

M^, = 2egpDDg'"'s2,, 
M%a = y/2egpDD{ki-2k3Yr^ 



t — mjj 



X (fci - ki + kAY£2u, 

A^^(, = V2e5^DD[(-2fci-fc2)V 

+ (fci + 2fc2)''5"'' + (fci - fc2)"5n 



s — m„ 



9al3 



(fci + fc2)c<(fcl + fc2)/3 



X (fca - fc4)''ir2i., 



-\/2, 



egpDDg'^''£2i^, 



M^a = e5pD-i5-[(-fci-fc3)"5''^ + (2fci-fc3)V^ 

+ (2fc3-fcir5""]^4;^ 



t — mjj. 



Jaf3 



(fci - fc3)c<(fcl - fc3)/3 



X [-2p^.g^'^ + (fc2 + fc4)''<?"" - 2k^gf^] 
X e2i/e3A£4uj) 
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+ (fci-2fc4r5n- 



u — m. 



5a/3 



(fci - ki)aih - ki))} 



X [(-/e2-fc3)V^ + 2fc2V'^ + 2fc3V^] 
X £2L/£3A£4aj) 

M"^, = \/2e5,D.D.[(-fci-fc3r5''^ + (2pi-p3)V 

1 



+ (2/03 - fcifs"^] 



9oil3 



{ki - fc3)a(fci - fc3)/3 



X e2L/£3A£4aj, 

A^^^ = -V2egpD'D' [{-2ki - fe)"^''" + (fci + 2A;2)''ff" 

+ (fci-fc2r5n- 



— 



(fcl + k2)a{kl + k2)0 



X [-2/c^/'^ + 2fc^/^ + (A;4- A;3)V'^] 

X £2ve3\eAuj- (5) 

The cross sections for these reactions with three par- 
ticles in the final state can be expressed in terms of the 
off-shell cross sections of the subprocesses involving two 
particles in the final state. Following the method of Ref. 
[23] for the reaction NN — !■ NAK, the spin-averaged dif- 
ferential cross section for the four reactions in Fig. 1 can 
be written as 



da. 



-yN^ND'D 



dtds 



327r2sp2 



while those for the two reactions in Fig. 2 are 



(6) 



da. 



■yN^NDD 



dtds\ 



6477^ sp^ 



(t - m2^2 



[4(1 + K,)2 



2 +\2 



+ 4(l-f Kp) 



X Kpi'^ml - t)] ap^^jjc{s^,t)\F{t)\''. (7) 

In the above, Pc is the ccntcr-of-mass momentum of 7 
and N, t is the squared four momentum transfer, and 
s is the squared center-of-mass energy. The quantity Si 
and k arc, respectively, the squared invariant mass and 
center-of-mass momentum of the tt and 7 in the subpro- 
cess 7A'' — > D*DN or of the p and 7 in the subprocesses 



7iV DON and 7iV D*D*N. Cross sections for 
these subprocesses are obtained from the amplitudes in 
Eqs.(4) and (5) using the software package FORM [24] 
to evaluate the summation over the polarizations of both 
initial and final particles. The differential cross sections 
for the two reactions 7iV D*D*N in Fig. 3 are similar 
to those for jN DDN with ap,y^jjij{si,t) replaced by 
ap^^jy,Q,{si,t). 

We have introduced the form factors F^jv jv and Fpjvjv 
at the ttA'^A'' and pNN vertices, respectively, to take into 
account the finite size of hadrons. As in Ref. [9], both 
are taken to have the following monopole form: 



F{t) 



A2 



m 



(8) 



where m is the mass of the exchanged pion or rho meson, 
and A is a cutoff parameter. Following Ref. [9] we take 
AjrATAr = 1.3 GcV and Ap^N = 1.4 GeV. We have also in- 
troduced a universal form factor at the strong interaction 
vertices in the ttj — > D*D, DD, and pj D*D* 

two-body subprocesses. Such a prescription guarantees 
the gauge invariance when all diagrams are included in 
each process [25]. The form factor we include here is of 
the following dipole form: 



A2 



A2 



(9) 



with q denoting the three momentum of photon in the 
center-of-mass system. We choose the cutoff parameter 
A to be 1.45 GeV to best reproduce the data. This form 
is also our default choice of form factors in this paper, 
unless stated otherwise. 

The cross sections for charmed hadron production 
from neutron is the same as that from proton if we ne- 
glect diagrams involving photon coupling directly to pro- 
ton. The isospin-averaged cross section for the reaction 
jN DDN is thus given by the sum of the cross sec- 
tions for the four processes in Fig. 1, which are ob- 
tained by integrating Eq.(6) over t and si. Similarly, 
one can obtain from Eq.(7) the isospin averaged cross 
sections for the reactions 77V DD*N{D*DN) and 
jN — !■ D*D*N, shown respectively, in Figs. 2 and 3. 
The results are shown in Fig. 4 by the solid, dotted, 
and dashed curves, respectively, for the reactions 7iV — > 
DDN, jN ^ DD*N{D*DN) and 77V D*D*N. It 
is seen that the reaction 7iV DD*N{D*DN) has the 
largest cross section with a peak value of about 40 nb, 
while the reaction 7A'' DDN has the smallest cross 
section of only about 1 nb. The larger cross sections for 
processes involving a charmed vector meson in the final 
state is due to the presence of interaction vertices with 
three vector mesons, which have a stronger momentum 
dependence than vertices with fewer number of vector 
mesons, leading thus to a larger strength at high ener- 
gies. 
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FIG. 4. Cross sections for photoproduction of charmed 
hadrons from nucleon with three particles in the final states. 

The above results are obtained without contributions 
from diagrams involving the photon coupled directly to 
external nucleons. These diagrams are needed to pre- 
serve the gauge invariance in each process. Their contri- 
butions are small compared to those from diagrams with 
the photon coupled directly to mesons. This is due to 
the s-channcl nucleon propagator (l/(s — n^'pf)) in these 
diagrams, which suppresses their amplitudes more than 
the t-channel heavy meson propagator in other as a result 
of the large photon energy needed to produce both the 
charmed and anticharmed meson pair. In the following, 
we demonstrate this effect by comparing the contribution 
due to diagram (Id) with that due to diagrams (la)-(lc) 
in Fig. 1. 

The amplitude for diagram (Id) in Fig. 1 can be writ- 
ten as 



M = i2eg„NNg7vDD' 



1 



{s - m%){t - ml) 

X P(P3)75(A + ^2 + '>^Nh^p{Pl)£nP5Su 
= 2g77DD'M2P^S„, 



(10) 



where pi, p3, p2, and are the momenta of initial and 
final nucleons, photon, and charmed meson, respectively. 
The cross section due to this diagram alone is given by 



da 



dtdsi 2567r^spg 



|A^2rr(si), 



(11) 



where again Pc is the center-of-mass momentum of the 

nucleon and the photon, si is the invariant mass of the 
£)*-£)+ pai]-, and r(si) is the decay width of the off-shell 
^0 ^ D*-D+. 
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FIG. 5. Cross sections due to photon coupling directly to 

nucleon (diagram (Id) in Fig. 1, dashed curve) and due to 
pion exchange (diagrams (la)-(lc) in Fig. 1, solid curve). No 
form factors are included at interaction vertices. 

The cross section due to the s— channel diagram (Id) 
in Fig. 1 involving photon coupling directly to nuclcion is 
shown by the dashed curve in Fig. 5 together with that 
coming from the pion exchange contributions (diagrams 
(la)-(lc) in Fig. 1), shown by the solid curve. Form fac- 
tors have been neglected in these results as we are only 
interested in their relative magnitude. It is seen that 
the contribution from the diagram with direct photon- 
nucleon coupling is more than two orders of magnitude 
smaller than that from the pion-exchange diagrams and 
can thus be safely neglected. 



III. PHOTOPRODUCTION OF CHARMED 
HADRONS FROM NUCLEON WITH 
TWO-BODY FINAL STATES 



Charmed hadron can also be produced in photo- 
nucleon reaction from processes involving two particles in 
the final state, i.e., — > DKc and D*Ac, shown 

by diagrams in Fig. 6. The interaction Lagrangians 
needed to evaluate the cross sections for these reactions 
are: 

Cd'na^ = gD^NAANli.KD*'' + D*^'K^^N). (12) 

As in ref. [15], we use SU(4) relations to determine the 
coupling constants goNAc and go'NAc in terms of known 
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coupling constants q^nn and QpNN, and they are given 
by 



3 - 2aD ^ „ . 

9DNAc = 7= QttNN — 9-kNN = io.o, 



V3 

-^f^gpNN = -5.6, 



(13) 



where ao = D/{D + F) ~ 0.64 [26] with D and F be- 
ing the coefficients for the usual D-type and F-type cou- 
plings. Values of these couplings constants are similar to 
those obtained from a QCD sum rule analysis [27]. 
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FIG. 6. Photoproduction of charmed hadrons from nucleon 
with two-particle final states. 

The amplitudes for the two processes in Fig. 6 are 



M^2 = {M'i^a + Mt^,)e^e^, 



(14) 



with A^9„, A^9f,, Mwa, and A^ioc f"^' ^^"^ ^'^^^ 
grams in Fig. 6, while M'^'^^, M^"^^, M'^^^, and M'^^,, for 
for the bottom four diagrams. They are given explicitly 

by 



■^9a = iegoNAc 



t — m|j 



{P2 - 2piYKc{pz)l^n{p^), 



iegDNA^ 



u — rriK 



i '"-TV 

X Ac(p3)7^(A + ^ + TOjv)7''p(Pi), 
■^106 = iegDNA 



u — mt 



Mtia = egD'NAAc{P3h°'n{pi)- — 



-9al3 + 



{Pl-P3)a{Pl-P3)0 



m 



D* 



1 



■^m = ego'NA 



■ u — mi; 



X Ac(p3)7''(A -h + mj^ch''n{pi), 

■^12a = e.to-AfA, T 

s — mj^ 

X K{P3)Y{jt)i +i>2+^N)l'"p{pi), 

■^126 = ego-NA,- — -2- 

X Ac(P3)7''(A - 4 + "IaJ7''p(Pi)- 



(15) 



Here, pi, p2, Pa, and p4 denote the momentum of 7, A'', 
D{D*), and Ac, respectively. 

The spin-averaged differential cross sections for the 
four processes in Fig. 6 are then 



da. 



dt 



1 



2567rsp; 



:\M.\'\m'r,, 



(16) 



where i = 9, 10, 11, 12 for the the four reactions 
D-A+, jp 50A+, jn D*-h + , and 7 ^ D-K+. 

The isospin-averaged cross sections for the two reac- 
tions DAc and DAc can be obtained from 
the above cross sections, and they are given by 
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FIG. 7. cross section photoproduction of charmed hadrons 
via charmed meson exchange. 

Taking the form factor at interaction vertices to be of 
the dipole type as in Eq. (9) and with a cutoff param- 
eter A = 1.45 GcV, which is chosen to reproduce the 
experimentally observed relative strength between two- 
body and three-body decay final states in photoproduc- 
tion of charmed hadrons on nucleon [16], the cross sec- 
tions for the reactions 7A'' DAc (dashed curve) and 
'jN D*Ac (solid curve) are shown in Fig. 7 as func- 
tions of total center-of-mass energy. It is seen that both 
cross sections are large with a peak value of about 30 nb 
for 7iV DAc and about 50 nb for 7 A/' ^ D^A^. 



IV. TOTAL CROSS SECTION FOR CHARMED 
HADRON PRODUCTION IN 
PHOTON-NUCLEON REACTION 



J3 




s^'' (GeV) 

FIG. 8. Total and partial cross sections for charmed hadron 
production in photon-nucleon reactions as functions of cen- 
ter-of-mass energy. 

The total cross section for photoproduction of charmed 
hadrons on nucleon is given by the sum of the cross sec- 
tions for two-body and three-body final states. In Fig. 8, 
we show the total cross section (solid curve) together with 
that for two-body final state (dotted curve) and three- 
body final state (dashed curve). It is seen that two-body 
final states involving Ac and a charmed meson dominates 
at low energy, while the three-body final state involving 
a nucleon as well as a charmed and anticharmed meson 
pair is more important at high energies. The two have 



comparable magnitudes around center-of-mass energy of 
about 5.7 GeV. This is consistent with the experimental 
data at 6 GeV [16], which shows that the final state with 
a Ac constitutes about 35% of the total cross section. 

Charm production from photo-nucleon reaction can 
be estimated using the leading-order perturbative QCD 
[28-30]. The cross section in this approach is given by 



a 



-yN 



J2n 



2ml/ V 



dx (j'^^{vx)g{x), 



(18) 



where rric is the charm quark mass, g{x) is the gluon 
distribution function inside the nucleon, and v ~ p ■ 
with p and p-^ being the momenta of the incoming nu- 
cleon and photon. The cross section a'^3{(jj) is that for 
charm- anticharm quark production from the leading or- 
der photon-gluon scattering, i.e., 



{UJ) 



X log- 



9 a;2 

l + Jl 



1 + 



4m^ 



I + /1 

Ami 



- 1 + 



4m; 



Auril 



(19) 



where co^ = 2pg ■ p^, with the gluon momentum denoted 
by Pa- 
using rric = 1.3 GeV and the leading order MRST 2001 
parameterization of the gluon distribution function in nu- 
cleon [31], we have calculated the photo charm produc- 
tion cross section on nucleon using the LO QCD formula, 
and the result is shown Fig. 9 by the dashed curve to- 
gether with that from the effective hadronic model (solid 
curve) and the available experimental data (open circles). 
We see that the LO QCD result reproduces the data at 
6 GeV and at higher energies. However, the QCD pre- 
diction below 6 GeV falls well below the results from the 
effective hadronic model. It is known that the QCD for- 
mula for photoproduction of heavy quarks should work 
best when the momenta involved are of order the heavy 
quark mass rric. Below this momentum and near the 
threshold energy, large logarithms will appear in the per- 
turbative QCD approach and spoil its convergence [32]. 
At these energies, our phcnomenological hadronic ap- 
proach should be more reliable as the cross section is 
dominated by two-body final states with no additional 
contribution to cause any large correction. On the other 
hand, the results from the hadronic model at higher ener- 
gies fall short of the experimental data. This is expected 
because at these energies, contributions from four-body 
final state and from the exchange of heavier mesons will 
become important. At these higher energies, perturba- 
tive QCD calculations should be a more efficient way for 
determining the cross section for photo charm production 
than adding new ingredients into our phenomenological 
hadronic model. 
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shown by open circles. 



V. SUMMARY AND DISCUSSIONS 



In summary, the total cross section for charmed hadron 
production in photo-nucleon reaction is evaluated in an 
effective hadronic model. This model is based on a 
gauged SU(4) flavor symmetric Lagrangian with the pho- 
ton introduced as a C/cm(l) gauged particle. The symme- 
try breaking effect is taken into account via using empir- 
ical hadron masses and coupling constants. Form factors 
of the monopole type are introduced at interaction ver- 
tices. This hadronic model has been previously used to 
evaluate the dissociation cross section of J/ip by hadrons. 
For photo production of charmed hadrons on nucleon, we 
have included both two-body final states involving a A^. 
and a charmed meson as well as three-body final states in- 
volving a nucleon and a charmed and anticharmed meson 
pair. It is found that reactions with two-body final states 
dominate the production cross section at low center-of- 
mass energies while reactions with three-body states are 
more important at high center-of-mass energies. Using 
cut-off parameters in the form factors from previous stud- 
ies of J/tp absorption by hadrons, the model reproduces 
the lowest available experimental data at centcr-of-mass 
energy of 6 GeV. Our results thus again confirm the va- 
lidity of the effective hadron model in previous studies 
and provide a useful model for further studies of reac- 
tions involving heavy quarks at low and near threshold 
energies. 
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Charm production from photo-nucleon reaction in a hadronic model 
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We study the total cross section for photo production of charmed hadrons near threshold using 
a hadronic Lagrangian. Both two-body final states involving Ac and a charmed meson as well as 
three-body final states involving nucleon and a charm-anticharm meson pair are included. With 
appropriate cut-off parameters in the form factors at interaction vertices, the model gives a total 
cross section that is consistent with the measured data at center-of-mass energy of 6 GeV. The result 
is compared with the prediction from the leading-order perturbative QCD. 
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I. INTRODUCTION 



Reliable estimates of the production and scattering 
cross sections of open and/or hidden charmed hadrons 
in hadronic matter are important for understanding 
many phenomena in relativistic heavy ion collisions and 
hadron-nucleus reactions. In particular, the dissociation 
cross sections of J/4' by light meson or nucleon are di- 
rectly related to the role of hadronic suppression of J/^p 
production in heavy ion collisions [1,2], and have thus 
been a subject under active investigations. In the heavy 
charm quark mass limit, a formula based on the leading- 
order (LO) perturbative QCD was derived in Refs. [3,4] 
for the absorption cross sections of J/ip by hadrons. 
Higher-order corrections due to the target mass and rel- 
ativistic effects have also been estimated [5-7]. However, 
the perturbative QCD approach is useful only at very 
high energy when nontrivial higher twist and higher 
corrections are small. At low energies, these corrections 
become large, so nonperturbative approaches are needed. 
Various phenomenological models have thus been intro- 
duced for studying the J/4' absorption cross section by 
hadrons at low energies. These include the hadronic 
model based on effective Lagrangians [8-10], the QCD 
sum rules [11], and the quark-exchange model [12]. Re- 
sults from these phenomenological models all give much 
larger cross sections for J/tp absorption by light mesons 
than that given by the perturbative QCD formula. To 
test the prediction from the hadronic model, the same 
effective Lagrangian has been used to evaluate the ab- 
sorption cross section of J/ip by nucleon, and it is found 
to be consistent with that extracted from photonuclear 
production of J/ ip [9] . 

In the present paper, we generalize the effective 
hadronic Lagrangian, which has also been used to 
study charmed meson scattering by hadrons [13,14] and 
charmed meson production from meson-nucleon scatter- 
ing [15], to include the photon and to study charmed 
hadron production from photo-nucleon reaction near 
threshold. Both two-body {AcD) and three-body (NDD) 
final states are included. We find that using reason- 
able values for the cutoff parameters at the interaction 
vertices, the resulting charmed hadron production cross 



section near threshold energy is consistent with the mea- 
sured one at center-of-mass energy of 6 GeV [16], al- 
though at high energies its value is much smaller than 
that measured experimentally or given by the LO per- 
turbative QCD. We further find that the relative con- 
tribution of two-body to three-body final states is also 
consistent with the experimental data. Our results thus 
provide another independent test and confirmation of the 
validity of hadronic model for determining the cross sec- 
tions involving heavy flavors at low energies. The ef- 
fective hadronic Lagrangian used in the present study 
will also be useful for evaluating other low energy cross 
sections involving heavy flavored hadrons, which can be 
studied at both the Japanese Hadron facility [17] and at 
the GSI future accelerator [18], where open and hidden 
charmed hadrons are copiously produced in proton- and 
antiproton- nucleus reactions near threshold. 

This paper is organized as follows. In Section II, we 
first calculate the cross sections for photoproduction of 
charmed hadrons from nucleon with three particles in the 
final state, which is dominated by pion and rho meson ex- 
change. Photoproduction of charmed hadrons from nu- 
cleon with two particle in the final state, which includes 
the charmed meson exchange, is studied in Section HI. 
In Section IV, we show the total cross section and com- 
pare it to available experimental data and results from 
the LO QCD. Finally, summary and discussions are given 
in Section V. 



II. PHOTOPRODUCTION OF CHARMED 
HADRONS FROM NUCLEON WITH 
THREE-PARTICLE FINAL STATES 



In photoproduction of charmed hadrons from nucleon 
with three particle in the final states, the three possible 
reactions are 77V DDN, jN DD*N{D*DN) ,and 
7A^ — > D*D*N . The lowest-order diagrams for the pro- 
cess DDN are shown in Fig. 1 and involves the 
exchange of pion in the intermediate state; those for the 
processes 77V DD*N{D*DN) and 7iV D*D*N in- 
volve the exchange of rho meson, and the lowest-order 
diagrams for the two processes are shown in 2 and 3. 
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FIG. 1. Photoproduction of charmed hadrons from nucleon 
involving pion exchange. 
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FIG. 2. Photoproduction of charmed hadrons {DD) from 
proton involving rho meson exchange. 



To evaluate the cross sections for these processes, we 
use the same Lagrangian introduced in Refs. [9,14,15] for 
studying charmed meson scattering by hadrons. This La- 
grangian is based on the gauged SU(4) flavor symmetry 
but with empirical masses. The coupling constants arc 
taken, if possible, from empirical information. Other- 
wise, the SU(4) relations are used to relate the unknown 
coupling constants to the known ones. Photon is then in- 
troduced in the Lagrangian via gauging its Uem(l) part 
as in ref. [19]. 
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FIG. 3. Photoproduction of charmed hadrons {D*D*) from 
proton involving rho meson exchange. 

The interaction Lagrangian densities that are relevant 
to the processes shown in Figs. 1, 2, and 3 are given as 
follows: 



JCwNN = -ig-KNNN^^fN ■ TT 
^pNN 



jO-nDD' = i9nDD'D*'^T ■ {Ddf^TT - d^Dlf) + H.C., 

jCpD'D' = igpD'D' [{d^D*-fDl - D*-fd^Dl) ■ ^ 
+ {D*-f-d^fr-d^D*''f-p,)D*>^ 



t^p-^DD 
l^p-yD*D* 



= ieA^[DQd^,D - {d^D)QD], 

= ie[A^'{d^D*''QDl - D^^Qd^Dl) 

+ {d^A-Dl - A''d^Dl)QD*^' 

+ D*^^Q{A^'^^Dl-^^,A-Dl], 

= -eg^DD'A>'{D*^{2TQ - Qt)D 

+ D{2QT-TQ)D;)-n, 

= egpooAi'DiTQ + Qt)D ■ p^, 

= egpD'D'{A''Dt{2TQ - Qr)D; 

+ A-'D;{2rQ~QT)D: 

-A^D*-'{2tQ-Qt)D:)-p^. 



(1) 



In the above, t arc Pauli spin matrices, and tt and p de- 
note the pion and rho meson isospin triplet, respectively, 
while D = {D^,D+) and D* = {D*°,D* + ) denote the 
pseudoscalar and vector charm meson doublets, respec- 
tively. The operator Q is the diagonal charge operator 
with diagonal elements equal to and -1. 

For coupling constants in the above interaction La- 
grangians, we use the following empirical values: g^NN = 
13.5 [20], gpNN = 3.25, Hp = 6.1 [21], g^oD- = 5.56 
[15,22], and 5^z3Z3 = 2.52 [14]. 

The diagrams in Figs. 1,2, and 3 can be separated to 
two types; one in which the photon is coupled to mesons 
such as the first three diagrams (denoted by (ia) to (ic) 
with i=l to 8), and the other in which the photon is 



2 



coupled directly to the incoming or outgoing charged nu- 
cleon. As shown later, contributions from the latter type 
are much smaller than those from the first type of di- 
agrams and are thus neglected in the following calcula- 
tions. We note that diagrams of the first type are similar 
to those for J/ip absorption on nucleon, which can be 
interpreted as absorption by the virtual pion and rho 
meson cloud of a nucleon. Here, they can be considered 
as charmed hadron production from the nucleon's virtual 
meson cloud. 

The amplitudes for the four processes in Fig. 1 are 
given by 



Mi = -iag^NNN{p3)-/5N{pi)-^ 

X {M,a+Mib + M,^), 



1 



(2) 



with i = 1 to 4, while the amplitudes for the four pro- 
cesses in Fig. 2, 3 can be written as 



Mj = agpNNN(j>3) 



ipi -P3)a 
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t-mj 



(3) 



with j = 5 to 8. In the above, pi and are four mo- 
menta of the initial and final nucleons, respectively. The 
coefficient a is, respectively, 1 and \/2 for neutral and 
charged pion coupling to nucleon. 

The three amplitudes Mia, Mn,, and M^c are for the 
subprocess ivf — > D*D in Fig. 1. It can be shown that 
they fulfill the chiral constraint [10], i.e., their sum van- 
ishes at the soft pion limit. Explicitly, they are given 

by 
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(4) 



where ki denotes the momentum of particle i of each 
subprocess, and and Si, are the polarization vector of 
7 and D* , respectively. We choose the convention that 
particles 1 and 2 represent initial-state particles while 
particles 3 and 4 represent final-state ones on the left 
and right sides of the diagrams. 

The amplitudes MJ^, Mj^, and Mj^ are those for the 
subprocesses pj DD and pj D*D* in Figs. 2 and 
3, and they are are written explicitly as: 
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The cross sections for these reactions with three par- 
ticles in the final state can be expressed in terms of the 
off-shell cross sections of the subprocesses involving two 
particles in the final state. Following the method of Ref. 
[23] for the reaction NN — !■ NAK, the spin-averaged dif- 
ferential cross section for the four reactions in Fig. 1 can 
be written as 



da. 



-yN^ND'D 



dtds 



327r2sp2 



while those for the two reactions in Fig. 2 are 



(6) 
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■yN^NDD 



dtds\ 



6477^ sp^ 



(t - m2^2 



[4(1 + K,)2 



2 +\2 



+ 4(l-f Kp) 



X Kpi'^ml - t)] ap^^jjc{s^,t)\F{t)\''. (7) 

In the above, Pc is the ccntcr-of-mass momentum of 7 
and N, t is the squared four momentum transfer, and 
s is the squared center-of-mass energy. The quantity Si 
and k arc, respectively, the squared invariant mass and 
center-of-mass momentum of the tt and 7 in the subpro- 
cess 7A'' — > D*DN or of the p and 7 in the subprocesses 



7iV DON and 7iV D*D*N. Cross sections for 
these subprocesses are obtained from the amplitudes in 
Eqs.(4) and (5) using the software package FORM [24] 
to evaluate the summation over the polarizations of both 
initial and final particles. The differential cross sections 
for the two reactions 7iV D*D*N in Fig. 3 are similar 
to those for jN DDN with ap,y^jjij{si,t) replaced by 
ap^^jy,Q,{si,t). 

We have introduced the form factors F^jv jv and Fpjvjv 
at the ttA'^A'' and pNN vertices, respectively, to take into 
account the finite size of hadrons. As in Ref. [9], both 
are taken to have the following monopole form: 



F{t) 



A2 



m 



(8) 



where m is the mass of the exchanged pion or rho meson, 
and A is a cutoff parameter. Following Ref. [9] we take 
AjrATAr = 1.3 GcV and Ap^N = 1.4 GeV. We have also in- 
troduced a universal form factor at the strong interaction 
vertices in the ttj — > D*D, DD, and pj D*D* 

two-body subprocesses. Such a prescription guarantees 
the gauge invariance when all diagrams are included in 
each process [25]. The form factor we include here is of 
the following dipole form: 



A2 



A2 



(9) 



with q denoting the three momentum of photon in the 
center-of-mass system. We choose the cutoff parameter 
A to be 1.45 GeV to best reproduce the data. This form 
is also our default choice of form factors in this paper, 
unless stated otherwise. 

The cross sections for charmed hadron production 
from neutron is the same as that from proton if we ne- 
glect diagrams involving photon coupling directly to pro- 
ton. The isospin-averaged cross section for the reaction 
jN DDN is thus given by the sum of the cross sec- 
tions for the four processes in Fig. 1, which are ob- 
tained by integrating Eq.(6) over t and si. Similarly, 
one can obtain from Eq.(7) the isospin averaged cross 
sections for the reactions 77V DD*N{D*DN) and 
jN — !■ D*D*N, shown respectively, in Figs. 2 and 3. 
The results are shown in Fig. 4 by the solid, dotted, 
and dashed curves, respectively, for the reactions 7iV — > 
DDN, jN ^ DD*N{D*DN) and 77V D*D*N. It 
is seen that the reaction 7iV DD*N{D*DN) has the 
largest cross section with a peak value of about 40 nb, 
while the reaction 7A'' DDN has the smallest cross 
section of only about 1 nb. The larger cross sections for 
processes involving a charmed vector meson in the final 
state is due to the presence of interaction vertices with 
three vector mesons, which have a stronger momentum 
dependence than vertices with fewer number of vector 
mesons, leading thus to a larger strength at high ener- 
gies. 
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FIG. 4. Cross sections for photoproduction of charmed 
hadrons from nucleon with three particles in the final states. 

The above results are obtained without contributions 
from diagrams involving the photon coupled directly to 
external nucleons. These diagrams are needed to pre- 
serve the gauge invariance in each process. Their contri- 
butions are small compared to those from diagrams with 
the photon coupled directly to mesons. This is due to 
the s-channcl nucleon propagator (l/(s — n^'pf)) in these 
diagrams, which suppresses their amplitudes more than 
the t-channel heavy meson propagator in other as a result 
of the large photon energy needed to produce both the 
charmed and anticharmed meson pair. In the following, 
we demonstrate this effect by comparing the contribution 
due to diagram (Id) with that due to diagrams (la)-(lc) 
in Fig. 1. 

The amplitude for diagram (Id) in Fig. 1 can be writ- 
ten as 



M = i2eg„NNg7vDD' 



1 



{s - m%){t - ml) 

X P(P3)75(A + ^2 + '>^Nh^p{Pl)£nP5Su 
= 2g77DD'M2P^S„, 



(10) 



where pi, p3, p2, and are the momenta of initial and 
final nucleons, photon, and charmed meson, respectively. 
The cross section due to this diagram alone is given by 



da 



dtdsi 2567r^spg 



|A^2rr(si), 



(11) 



where again Pc is the center-of-mass momentum of the 

nucleon and the photon, si is the invariant mass of the 
£)*-£)+ pai]-, and r(si) is the decay width of the off-shell 
^0 ^ D*-D+. 
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FIG. 5. Cross sections due to photon coupling directly to 

nucleon (diagram (Id) in Fig. 1, dashed curve) and due to 
pion exchange (diagrams (la)-(lc) in Fig. 1, solid curve). No 
form factors are included at interaction vertices. 

The cross section due to the s— channel diagram (Id) 
in Fig. 1 involving photon coupling directly to nuclcion is 
shown by the dashed curve in Fig. 5 together with that 
coming from the pion exchange contributions (diagrams 
(la)-(lc) in Fig. 1), shown by the solid curve. Form fac- 
tors have been neglected in these results as we are only 
interested in their relative magnitude. It is seen that 
the contribution from the diagram with direct photon- 
nucleon coupling is more than two orders of magnitude 
smaller than that from the pion-exchange diagrams and 
can thus be safely neglected. 



III. PHOTOPRODUCTION OF CHARMED 
HADRONS FROM NUCLEON WITH 
TWO-BODY FINAL STATES 



Charmed hadron can also be produced in photo- 
nucleon reaction from processes involving two particles in 
the final state, i.e., — > DKc and D*Ac, shown 

by diagrams in Fig. 6. The interaction Lagrangians 
needed to evaluate the cross sections for these reactions 
are: 

Cd'na^ = gD^NAANli.KD*'' + D*^'K^^N). (12) 

As in ref. [15], we use SU(4) relations to determine the 
coupling constants goNAc and go'NAc in terms of known 
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coupling constants q^nn and QpNN, and they are given 
by 



3 - 2aD ^ „ . 

9DNAc = 7= QttNN — 9-kNN = io.o, 



V3 

-^f^gpNN = -5.6, 



(13) 



where ao = D/{D + F) ~ 0.64 [26] with D and F be- 
ing the coefficients for the usual D-type and F-type cou- 
plings. Values of these couplings constants are similar to 
those obtained from a QCD sum rule analysis [27]. 





n Y n T p T P 1^ 
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FIG. 6. Photoproduction of charmed hadrons from nucleon 
with two-particle final states. 

The amplitudes for the two processes in Fig. 6 are 



M^2 = {M'i^a + Mt^,)e^e^, 



(14) 



with A^9„, A^9f,, Mwa, and A^ioc f"^' ^^"^ ^'^^^ 
grams in Fig. 6, while M'^'^^, M^"^^, M'^^^, and M'^^,, for 
for the bottom four diagrams. They are given explicitly 

by 



■^9a = iegoNAc 



t — m|j 



{P2 - 2piYKc{pz)l^n{p^), 



iegDNA^ 



u — rriK 



i '"-TV 

X Ac(p3)7^(A + ^ + TOjv)7''p(Pi), 
■^106 = iegDNA 



u — mt 



Mtia = egD'NAAc{P3h°'n{pi)- — 



-9al3 + 



{Pl-P3)a{Pl-P3)0 



m 



D* 



1 



■^m = ego'NA 



■ u — mi; 



X Ac(p3)7''(A -h + mj^ch''n{pi), 

■^12a = e.to-AfA, T 

s — mj^ 

X K{P3)Y{jt)i +i>2+^N)l'"p{pi), 

■^126 = ego-NA,- — -2- 

X Ac(P3)7''(A - 4 + "IaJ7''p(Pi)- 



(15) 



Here, pi, p2, Pa, and p4 denote the momentum of 7, A'', 
D{D*), and Ac, respectively. 

The spin-averaged differential cross sections for the 
four processes in Fig. 6 are then 



da. 



dt 



1 



2567rsp; 



:\M.\'\m'r,, 



(16) 



where i = 9, 10, 11, 12 for the the four reactions 
D-A+, jp 50A+, jn D*-h + , and 7 ^ D-K+. 

The isospin-averaged cross sections for the two reac- 
tions DAc and DAc can be obtained from 
the above cross sections, and they are given by 
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FIG. 7. cross section photoproduction of charmed hadrons 
via charmed meson exchange. 

Taking the form factor at interaction vertices to be of 
the dipole type as in Eq. (9) and with a cutoff param- 
eter A = 1.45 GcV, which is chosen to reproduce the 
experimentally observed relative strength between two- 
body and three-body decay final states in photoproduc- 
tion of charmed hadrons on nucleon [16], the cross sec- 
tions for the reactions 7A'' DAc (dashed curve) and 
'jN D*Ac (solid curve) are shown in Fig. 7 as func- 
tions of total center-of-mass energy. It is seen that both 
cross sections are large with a peak value of about 30 nb 
for 7iV DAc and about 50 nb for 7 A/' ^ D^A^. 



IV. TOTAL CROSS SECTION FOR CHARMED 
HADRON PRODUCTION IN 
PHOTON-NUCLEON REACTION 



J3 




s^'' (GeV) 

FIG. 8. Total and partial cross sections for charmed hadron 
production in photon-nucleon reactions as functions of cen- 
ter-of-mass energy. 

The total cross section for photoproduction of charmed 
hadrons on nucleon is given by the sum of the cross sec- 
tions for two-body and three-body final states. In Fig. 8, 
we show the total cross section (solid curve) together with 
that for two-body final state (dotted curve) and three- 
body final state (dashed curve). It is seen that two-body 
final states involving Ac and a charmed meson dominates 
at low energy, while the three-body final state involving 
a nucleon as well as a charmed and anticharmed meson 
pair is more important at high energies. The two have 



comparable magnitudes around center-of-mass energy of 
about 5.7 GeV. This is consistent with the experimental 
data at 6 GeV [16], which shows that the final state with 
a Ac constitutes about 35% of the total cross section. 

Charm production from photo-nucleon reaction can 
be estimated using the leading-order perturbative QCD 
[28-30]. The cross section in this approach is given by 



a 



-yN 



J2n 



2ml/ V 



dx (j'^^{vx)g{x), 



(18) 



where rric is the charm quark mass, g{x) is the gluon 
distribution function inside the nucleon, and v ~ p ■ 
with p and p-^ being the momenta of the incoming nu- 
cleon and photon. The cross section a'^3{(jj) is that for 
charm- anticharm quark production from the leading or- 
der photon-gluon scattering, i.e., 



{UJ) 



X log- 



9 a;2 

l + Jl 



1 + 



4m^ 



I + /1 

Ami 



- 1 + 



4m; 



Auril 



(19) 



where co^ = 2pg ■ p^, with the gluon momentum denoted 
by Pa- 
using rric = 1.3 GeV and the leading order MRST 2001 
parameterization of the gluon distribution function in nu- 
cleon [31], we have calculated the photo charm produc- 
tion cross section on nucleon using the LO QCD formula, 
and the result is shown Fig. 9 by the dashed curve to- 
gether with that from the effective hadronic model (solid 
curve) and the available experimental data (open circles). 
We see that the LO QCD result reproduces the data at 
6 GeV and at higher energies. However, the QCD pre- 
diction below 6 GeV falls well below the results from the 
effective hadronic model. It is known that the QCD for- 
mula for photoproduction of heavy quarks should work 
best when the momenta involved are of order the heavy 
quark mass rric. Below this momentum and near the 
threshold energy, large logarithms will appear in the per- 
turbative QCD approach and spoil its convergence [32]. 
At these energies, our phcnomenological hadronic ap- 
proach should be more reliable as the cross section is 
dominated by two-body final states with no additional 
contribution to cause any large correction. On the other 
hand, the results from the hadronic model at higher ener- 
gies fall short of the experimental data. This is expected 
because at these energies, contributions from four-body 
final state and from the exchange of heavier mesons will 
become important. At these higher energies, perturba- 
tive QCD calculations should be a more efficient way for 
determining the cross section for photo charm production 
than adding new ingredients into our phenomenological 
hadronic model. 
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FIG. 9. Cross sections for charm production in photo nu- 

clcon reaction in the hadronic model (solid curve) and the 
pQCD approach (dashed curve). The experimental data are 
shown by open circles. 



V. SUMMARY AND DISCUSSIONS 



In summary, the total cross section for charmed hadron 
production in photo-nucleon reaction is evaluated in an 
effective hadronic model. This model is based on a 
gauged SU(4) flavor symmetric Lagrangian with the pho- 
ton introduced as a C/cm(l) gauged particle. The symme- 
try breaking effect is taken into account via using empir- 
ical hadron masses and coupling constants. Form factors 
of the monopole type are introduced at interaction ver- 
tices. This hadronic model has been previously used to 
evaluate the dissociation cross section of J/ip by hadrons. 
For photo production of charmed hadrons on nucleon, we 
have included both two-body final states involving a A^. 
and a charmed meson as well as three-body final states in- 
volving a nucleon and a charmed and anticharmed meson 
pair. It is found that reactions with two-body final states 
dominate the production cross section at low center-of- 
mass energies while reactions with three-body states are 
more important at high center-of-mass energies. Using 
cut-off parameters in the form factors from previous stud- 
ies of J/tp absorption by hadrons, the model reproduces 
the lowest available experimental data at centcr-of-mass 
energy of 6 GeV. Our results thus again confirm the va- 
lidity of the effective hadron model in previous studies 
and provide a useful model for further studies of reac- 
tions involving heavy quarks at low and near threshold 
energies. 
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